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ABSTRACT 



We have carried out a Spitzer survey for brown dwarf disks in the ~5Myr 
old Upper Scorpius (UpSco) star forming region, using IRS spectroscopy from 8 
to 12 fim and MIPS photometry at 24 fim. Our sample consists of 35 confirmed 
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very low mass members of UpSco. Thirteen objects in this sample show clear 
excess flux at 24 fim, explained by dust emission from a circum-sub-stellar disk. 
The spectral energy distributions (SEDs) of the remaining objects are consistent 
with pure photospheric emission. Objects without excess emission either have no 
disks at all or disks with inner opacity holes of at least ~ 5AU radii. Our disk 
frequency of 37 ± 9% is higher than what has been derived previously for K0-M5 
stars in the same region (on a 1.8a confidence level), suggesting a mass-dependent 
disk lifetime in UpSco. The clear distinction between objects with and without 
disks as well as the lack of transition objects shows that disk dissipation inside 
5 AU occurs rapidly, probably on timescales of < 10 5 years. For the objects with 
disks, most SEDs are uniformly flat with flux levels of a few mJy, well modeled 
as emission from dusty disks affected by dust settling to the midplane, which 
also provides indirect evidence for grain growth. The silicate feature around 
10 [im is either absent or weak in our SEDs, arguing for a lack of hot, small dust 
grains. Compared with younger objects in Taurus, brown dwarf disks in UpSco 
show less flaring. By comparing SEDs of stars and brown dwarfs in UpSco, we 
find that dust settling is not a strong function of mass in this region. Taken 
together, these results clearly demonstrate that we see disks in an advanced 
evolutionary state: Dust settling and grain growth are ubiquituous in circum- 
sub-stellar disks at ages of 5 Myr, arguing for planet forming processes in brown 
dwarf disks. For almost all our targets, results from high-resolution spectroscopy 
and high-spatial resolution imaging have been published before, thus providing 
a large sample of brown dwarfs for which information about disks, accretion, 
and binarity is available. We find that nine out of 13 objects with disks do 
not accrete significantly. Hence, dusty disks can persist although the continuous 
accretion has stopped or dropped below measurable levels. Four objects with 
disks are binaries, three of them with (projected) separations < 30 AU. These 
objects likely harbour small disks truncated by the companion. 

Subject headings: stars: circumstellar matter, formation, low-mass, brown dwarfs 
- infrared: stars 



1. Introduction 



Young brown dwarfs are known to undergo a T Tauri phase similar to solar-mass stars. 
A la rge fraction of them exhibits spectroscopic signatures of ongoing accreti on and outflows 



(.e.g. 



Jayawardhana et al.ll2003bl : iMohanty et al.ll2005l ; iMuzerolle et al.ll2005l ) combined with 
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strong, T Tauri like variability (e.g. IScholz fc Eisloffell 120051 ; IScholz &: Jayawardhanal 120061 ) . 
Moreover, near/mid infrared colour excess has been found, indicative of the presence of 
dusty circum-sub-stellar disks. It is this latter aspect of the T Tauri phase in brown dwarf 
evolution that we want to investigate in this paper. 

Probing the properties and evolution of brown dwarf disks by analyzing their spectral 
energy distribution (SED) in the infrared and sub/mm regime is a useful test of brown 
dwarf origins, for example to distinguish between in situ formation and formation as stellar 
embryos ejected from a multiple system by searching for signatures of disk truncation (see 
Whitworth et al.ll2006l . for a summary). Additionally, disks are believed to be the nurseries of 
planets, and hence substellar disks allow us to study the early stages of planet formation, i.e. 
the process of grain and planetesimal growth, in an extreme case, constraining the efficiency 
and universality of current scenarios for the formation of planetary systems. 

Until recently, our knowledge of substellar disks relied on a few case studies. In th e 
pre-Spitzer era, m i d-infrared SEDs have b een published for example by iNatta et al.l (120021 ); 
Apai et al.l (120021 ); iMohanty et al.l (j2004bl ). in total for ~ 15 bro wn dwarfs wi t h age s of 1- 
2Myr, but only for one older object (the ~ 8Myr old 2M1207, Isterzik et all (bo04h ). By 
comparing the observed spectral energy distributions (SEDs) with model predictions, it has 
been found that young brown dwarfs appear to harbour disks with a diversity of geometries, 
which are interpreted as a consequence of disk evolution, consistent with findings for T Tauri 
stars. In some cases the opening angle of the disks increase with increasing distance from the 
star ("flared disks"), in others the disk is flat, which is expla ined as a consequence of grain 



grow th and settling of the larger grains to the disk midplane ( lApai et al.ll2004j ; ISterzik et al. 
20041 ) . 



Larger samples have been investigated in the near-infrared, where ground -based ob- 
servations are much easier (e.g. iMuench et al.ll200ll ; iJayawardhana et al.ll2003al ). providing 
constraints on disk frequency as a function of age: Based on near- infrared (K-, L-band) 
colour excess, disk frequencies decline strongly between 1 and 5 Myr, and only very few 
brown dwarfs retain their disks longer than that. Thus, the lifetimes for brown dwarf disks 
are in the range of 5 Myr, which is comparable to the lifetimes of disks of solar-mass stars 
dHaisch et aDboOlh . 



It is known, however, that the near-infrared 
ies of brown dwarf disks 



le-g. 



4 urn) has major limitations for stud- 



Comeron et al.l 120001 ) . Problems in this wavelength regime 



include, for example, the low contrast between disk and photosphere, the restriction to the 
innermost parts of the disks (< 0.1 AU), and the presence of strong photospheric molecular 
bands for these cool objects. As a consequence, the disk census derived from near-infrared 
observations may be incomplete, and it is difficult to obtain reliable constraints on the disk 
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lifetime. Mid-infrared observations give much better sensitivity for disk detections, because 
the photospheric flux rapidly declines at A > 4 pm. The wavelength range 4 — 25 pm is ideal 
to probe flaring geometry and inner disk hole sizes by comparing SED with model predic- 
tions. Moreover, an analysis of shape and intensity of the silicat e feature around 1 pm is 
able to give detailed constraints on the warm disk chemistry (e.g. iMeeus et al.ll2003l ). 



Spitzer provides unprecedented sensitivity between 3 and 24 pm, and thus allows us 
for the first time to carry out effi cient mid-infrared observations of large object samples. 
Using Spitzer's IRS spectrograph, lApai et al.l (120051 ) find silicate emission in six substellar 
disks at ages of 1-2 Myr, and confirm previous claims for grain growth and crystallization. 
Spitzer also allows to push the limits i n disk detection in terms of object age and mass (e.g. 
Riaz et alJliooa iLuhman et aDl2005ah . 



Only broad wavelength coverage and large object samples can overcome the large de- 
gree of uncertainties when deriving disk properties, considering the numerous degeneracies 
inherent in radiation transfer modeling of dusty disks, and thus provide reliable constraints 
on the disk parameters. This consideration is the main guideline for our long-term program 
to characterize brown dwarf disks. In a first paper, we have combined 1.3 mm measurements 
for 20 young brown dwarfs in Taurus with Spitzer IRAC/MIPS photometry, to co nstrain 
fundamental parameters of substellar dis ks in a systematic way (jScholz et al.l 120061 ) . Our 
results, together with previous findings by iKlein et al.l (120031 ). show that brown dwarf disks 
can have masses of several Jupiter masses; at least 25% have outer radii > 10 AU. In Taurus 
there is no evidence for truncated disks due to an ejection process early in the life history of 
brown dwarfs, implying that most sub-stellar objects probably form in isolation. 

In this second paper, we focus on evolutionary effects like disk dissipation, inner disk 
clearing, dust settling, and grain growth, which are believed to be consequences or prerequi- 
sites of planet formation. The occurence of these processes in brown dwarfs disks has already 
been confirmed, mostly at ages of 1-2 Myr. The goal of this paper is to probe disk evolution 
for a large sample of somewhat older objects. The ideal target region for this project is the 
star forming ass ociation Upper Scorpius (Up Sco): With an age of 5 Myr and supposedly low 



age spread (e.g. iPreibisch fc Zinneckerl 119991 ). it is significantly older than the well-studied 
star forming regions in Chamaeleon, p Ophiuchus, and Taurus, the main targets in previous 
brown dwarf disk surveys, and thus presents disks at an evolved stage. Moreover, UpSco is 
nearby (distance ~ 145 pc) and < 



(Ardila et ali 


2000; 


Martin et al. 


2004) 



(< 10) of brown dwarfs in UpSc o have already bee n studied in the near-/mid infrared by 
Jayawardhana et al.l (j2003al ) and iBouy et al. (l2006ah . 



Using Spitzer IRS spectroscopy and MIPS 24 pm photometry, we carried out a disk sur- 
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vey for a large sample of 35 UpSco brown dwarfs. This dataset has been complemented by 
Ha measurements and results from companions surveys (mostly from the literature), which 
gives us the opportunity to probe connections between disks, accretion, and binarity. In Sect. 
[2] we introduce the sample, describe the observations, the data reduction, and the extraction 
of SEDs. Subsequently, we give a general discussion of our SEDs, and derive an estimate of 
the disk frequency based on MIPS photometry (Sect. [3]). In the subsequent section, we com- 
pare the observed fluxes with models to constrain disk geometry (Sect. 14. ip and chemistry 
(Sect. 14. 2p . The interpretation of these results is given in Sect. \5\, including a discussion of 
the accretion/disk and disk/binarity connection. Finally, we give our conclusions in Sect. [6j 



2. Targets, observations, data reduction 



This paper is based on Spitzer observations of 35 very low mass objects in the Upper 
Sc orpius (UpSco) sta r for ming region. The ta rgets have been selected from the surveys 



of Ardila et al. 



rl2oooh an d iMartm et al.l (120041 ). excluding objects later identified as non- 
members (iMohanty et al.ll2004ah . From these spectroscopically confirmed young members 



of UpSco, we select ed 35 objects with spectral type >M5, cor responding to effective tem- 
peratures 5,3000 K (IMohanty et al.l l2004ai iLuhman et al.l 120031 ) and thus masses < 0.1 M & 
( IBaraffe et al.lll998l ). Although some of these objects may have masses higher than the sub- 
stellar limit at 0.08 M@, for the sake of simplicity we will call them brown dwarfs in the 
following. For a more detailed assessment of object masses see Sect. 13.11 Our target sample 
thus comprises almost all brown dwarfs in UpSco known at the time of the observations 
and is listed in Table 13.11 It does no t include the very low mass objects in this region 
found recently by lSlesnick et al.l (120061 ) . We note that the objects in our sample are mostly 
unbiased with respect to their disk environment, since they have been identified based on 
colour/magnitude diagrams and spectroscopic youth signatures. 



A note on nomenclature: Objects from lArdila et al.l (120001 ) are called 'UScoCTIO XXX' 
in the discovery paper, where XXX is an identi fication numb e r; we will abbreviate them 
as uscoXXX in the following. The targets from IMartm et al.l (120041 ) have originally been 
detected based on DENIS photometry and are thus called 'DENIS-P J' plus sky coordinates 
in DENIS (e.g., DENIS-P J155556.0-204518.5). For brevity, we will name these objects 'usd' 
followed by the first six digits of the sky coordinates, e.g. usdl55556. To avoid ambiguities, 
we include the coordinates of all sources in Table [TJ 

All 35 objects were targets of Spitzer GO program #20435 (PI: Ray Jayawardhana) . 
Since six of these objects overlapped with the GTO program #248 (PI: Subhanjoy Mohanty), 
they have been observed in the framework of the GTO campaign. For all objects we obtained 
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IRS spectroscopy from 7.4 to 14.5/im (module SL, 1st order) and MIPS images at 24/im. We 
aimed for a sensitivity of ~0.2 mJy, providing clear detections of the photospheric continuum 
at 10/xm even for our faintest targets and excellent upper limits at 24/im, allowing us a 
reliable distinction between objects with and without disks. The objects were observed with 
MIPS in one standard dither cycle with 10 sec exposure time per position and small offsets, 
resulting in a total on-source time of ~ 160 sec. The only exceptions were the six targets 
observed in program #248, in these cases three cycles were carried out, in total ~ 480 sec 
on-source time. IRS observations for all objects were carried out with two cycles, two nod 
positions, and 60 sec ramp duration, resulting in a total on-source time of 240 sec. The peak- 
up accuracy (i.e. the difference between pointing and requested position) was typically 0'/2 
or less; most of this offset is along the slit and thus slit losses are negligible. 

We used the post-BCD pipeline products for further reduction (pipeline version S13.2.0). 
Photometry on the MIPS post-BCD images was carried out using daophot aperture photom- 
etry in IRAF with a large aperture of 13" and a sky annulus with inner and outer radii of 
20" and 32", respectively. After applying the aperture correction given in the MIPS data 
handbook, version 3.2.1, fluxes were converted to mJy. Fifteen brown dwarfs are clearly 
detected in the MIPS observations, with flux levels between 0.35 and 10.7 mJy. For the 
non-detections, we estimated a generic upper limit from our images by measuring fluxes in a 
large number of regions containing only sky and no objects. These regions were chosen to be 
close to the positions, where the objects are expected to be. After excluding some outliers, 
these measurements scatter between -0.2 and 0.2 mJy. This is the flux level obtained from 
pure noise, and we thus infer a 2a upper limit of 0.4 mJy for our MIPS non-detections. We 
note that this limit may vary somewhat from object to object due to varying background 
level, but we prefer to use this statistically robust estimate instead of deriving upper limits 
for each object individually based on small number statistics. We estimate the uncertainty 
of our MIPS fluxes to be ~ 10%, composed of uncertainties in calibration (~ 5%), aperture 
correction, and instrumental magnitudes. (According to the suggestions in the MIPS data 
handbook, colour corrections are negligible for our objects.) All 24/im fluxes are listed in 
TableU 

For the IRS spectra, we started with the coadded spectra from the BCD data products, 
pipeline version S13.2.0, and corrected for the sky background by subtracting the two nod 
positions. Extraction of the spectra was done using the Spitzer IRS Custom Extractor Spice, 
version 1.3. All our targets are visible as profile peak along the slit, albeit some of them at 
the detection limit; flux levels at 10 /im are between 0.2 and 8.5 mJy. Extraction was carried 
out with standard settings for the two nod positions, followed by averaging these two single 
datasets to produce the final spectrum for the respective object. Additionally, we extracted 
the sky noise at two positions in the slit not affected by the objects. After excluding outliers, 
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these noise spectra provide us with an uncertainty estimate for each datapoint in the science 
spectra. Typical uncertainties of the individual datapoints are in the range of 0.2-0.3 mJy, 
where 7 spectra are affected of excessive noise in the range of ~ 0.5 mJy, which can in 
most cases be attributed to background emission from the molecular cloud. Following the 
information in the IRS Data Handbook, version 2.0, we expect additional flux calibration 
uncertainties of ~ 10%. Because of low flux levels in most spectra and thus low signal-to- 
noise ratio, we binned all IRS spectra to a resolution of ~ 0.25 fim. This reduces the number 
of datapoints per spectrum from ~ 100 to ~ 20 and improves the signal-to-noise ratio by 
a factor of ~ 2. All analysis using the IRS data is done on these binned spectra (see Sect. 

Q. 

All IRS spectra were truncated to the wavelength range between 7.5 and 12.5 \im. Fluxes 
at longer wavelengths are increasingly affected by the '14 micron teardrop' effect, which is 
believed to be a type of scatter light, as described in the IRS data handbook. We calculated 
broadband fluxes at 9, 10, and 11 fim from the IRS spectra by averaging over a wavelength 
range of 1/im. Table [1] contains these values together with the average noise level for the 
individual datapoints in the IRS spectrum for all objects. 

Final spectral energy distributions (SEDs) were produced by combining IRS and MIPS 



Ardila et al. 


(2000 


) and 


Martin et al. 


(2004) 



tern (jCarpenterl l200ll ) . magnitudes were converted to fl uxes using the zeropo ints given by 



Skinner (1996), which are based on the CIT zeropoints in lBeckwith et al.l (119761 ). Vice versa, 
we calculated K — 24 \im colours by converting the MIPS fluxes to magnitudes using the 
zeropoint of 7.14 ± 0.0815 Jy given in the MIPS data handbook. In Fig. [I] we show all SEDs 
for our sources, normalized to the J-band fluxes. 

The brown dwarf usdl60603 appears to have a faint neighbour in the MIPS image, at a 
separation of 10" . Since this object is not seen in 2MASS images, we assume that this is a 
spurious source, most likely an asteroid. Because the separation is considerably larger than 
the slitwidth of IRS (3.7"), the IRS spectrum of usdl60603 is probably not affected by this 
neighbour. (Moreover, the spectrum has been taken three weeks before the MIPS image, 
thus it is unlikely that the neighbour was at the same position, if it is indeed an asteroid.) 
The large aperture of the MIPS photometry, however, does contain both objects. To correct 
for this contamination, we measured relative fluxes of both objects using a small aperture, 
and corrected the total flux accordingly. We obtain a flux of 2.84 mJy for the brown dwarf 
and 1.50 mJy for the (likely asteroidal) neighbour. 

We aimed to include Ha measurements for our targets in this paper. Of particular 



interest is the Ha 10% width, which is correlated with the disk accretion rate (INatta et al. 
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2004 ) and can therefore be used as a signature for accre t ion. For 19 of our obje c ts, H a 



10% widths have been published by iJayawardhana et al.l (120021 ) or iMohanty et al. J2005h . 
We observed the remaining 16 objects using the high-resolution spectrographs UVES at 
ESO/VLT as part of the program 077.C-0323(B) and MIKE at the 6.5 m Clay Magellan 
telescope on Las Campanas. For UVES we used dichroic 2 with central wavelengths of 437 
and 760 nm and a 0'/8 slit, resulting in R ~ 50000. For MIKE we used a slit width of 0.7", 
2x2 binning, and achieved a typical resolution of R ~ 25000. Data reduction and line 
measurements were carried out in a standard fashion, see IScholz fc Jayawardhanal (120061 ) for 
details. All Ha 10% widths are listed in Table [TJ Most of our targets have b een observed 



recen tly with high spatial resol ution to search for companions, either with HST (IKraus et al. 



20051 ) or with adaptive optics (IBouy et al.ll2006bl ). We added the results from these surveys 
to Table OQ 



3. Spectral energy distributions 

The aim of this section is to provide an overview of the SEDs for our targets and to 
determine for how many objects we see evidence for excess mid-infrared emission due to 
disks. It turns out that based on a comparison of the Spitzer 24/zm fluxes with photospheric 
fluxes expected for our targets, the disk frequency in our sample is 37±9%. We will start with 
deriving fundamental parameters for the objects, to be able to estimate photospheric fluxes. 
Subsequently, we will discuss the SEDs and probe if disk excess is required to interpret the 
mid-infrared SEDs. 



3.1. Fundamental parameters for the targets 



A reliable distance of Upper Scorpius has been determined based on the Hipparcos 
measurements for the high-mass stars in this association. On average the distance is 145 pc 
( Ide Zeeuw et al.lll999l). with a spread of ±15 pc, consistent with a roughly spherical shape 
for the association (IPreibisch et al.ll2002l ). 



The age of the pre-main sequence population in Upper Scorpius has been subject of 
several recent studies. Since we are interested in evolutionary effects, a reliable estimate of 
the age of UpSco is relevant. For the high-mass stars in this region, an age of 5-6 Myr has 



been found from several different m ethods, with no evidence for age spread (Ide Geus et al. 



19891 ). For the low- mass population, IPreibisch &: Zinneckerl (119991 ) estimate an age of 5 Myr 
based on an HR diagram, in agreement with the results for high-mass stars. According to 
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their analysis, the age s pread is not si g nifica nt and probably smaller than 2 Myr. This view 
has been challenged by ISartori et al.l (120031 ). who find an average age of 8-10 Myr with a 
large spread of 1-20 Myr for Upper Scorpius. Their HR diagram analysis, however, is mainly 
based on early-type stars, for which 1-20 Myr isochrones quickly converge to the ZAMS and 
thus do not allow a clear discrimination between different ages. Therefore, we believe they 
overestimate the age spread. 



Recently, iMohanty et al.l (j2004al ) measured gravities and temperatures from high-resolution 
spectroscopy for 11 of our targets. When compared with evolutionary tracks, 8 objects are 
clustered arou nd the 5 Myr isochrone, but three of them (uscol04, 128, 130) appear to have 
ages < 1 Myr. IMohanty et al.l (j2004al ) discuss in detail several scenarios to explain this dis- 
crepancy, and arrive at the conclusion that a large age spread is unlikely. Instead, they favour 
a scenario where this apparent age offset arises from evolutionary model uncert ainties related 



to acc retion, deuterium burning and/or convection effects. Thus, except for ISartori et al. 
(120031 ) , all recent studies give an age of ~ 5 Myr and no convincing evidence for strong age 



spread. This is the assumption which we will use in the remainder of this paper. Still there 
are some concerns, and a future study of ages in UpSco based on Lithium absorption is 
desirable. 



Visual e x tincti ons have been determined for our sample by lArdila et al.l (120001 ) and 
Martin et al. (j2004 ). They find an upper limit of Ay ~ 2 mag, for most objects the optical 
extinction is < 1 mag. This shows that a) near-infrared magnitudes are not significantly 
affected by extinction, and b) there is no indication for close to edge-on disks from the 
optical magnitudes. Therefore, we can neglect the effects of disk inclination for most of the 
following analysis. 

In Table [Q, we list effective temperatures for our 36 targets. These values have been 
de termined by comparin g the spectral types wit h the effective temperature scale provided 
by iLuhman et al.l (120031 ). IMohanty et al.l (l2004al ) derived fundamental parameters for some 
of our objects by spectral fitting, and their effective temperatures tend to be lower than our 
estimates by as much as 200 K. This discrepancy reflects the difficulties in assigning a reliable 
effective temperature scale for these cool objects, but might also be affected by problems 
in spectral typing. The absolute uncertainties of the effective temperatures in Table [T] are 
thus probably in the range of ±200 K. Since they have been derived in a consistent way, the 
relative accuracy is certainly better. 

The masses for our targets can be constrained by c omparing the effective tempera- 
tures with the 5 Myr isochrones from iBaraffe et al.l (120031 ). With this approa c h, we o btain 
masses between ~ 0.02 and 0.12 M & . The masses derived by IMohanty et al.l (j2004al ) span 
a larger range, from 0.01 to 0.26 M , where the discrepancies are probably again due to 
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some problems of the evolutionary tracks at young ages (see lBaraffe et al.ll2002l ). Still, these 
estimates confirm that all our objects have very low masses, although some of them may not 
be substellar. 



3.2. SEDs: General remarks 



The mid-infrared SEDs for the sample of 35 brown dwarfs in UpSco are plotted in Fig. 
[TJ In this Section, we give an overview of the derived SEDs purely based on the observational 
results, without invoking models for comparison. 

From the appearance of the SED, we can clearly divide the objects in two classes: There 
are thirteen objects with clear excess at 24 jum, which also have the highest fluxes between 
9 and 11 fim. The remaining ones are either not detected or have very low fluxes in MIPS, 
and are in most cases faint in IRS. In these two types of SED we likely see objects with and 
without disks, roughly speaking. (We will show in Sect. 13.31 more definitively that the 13 
with the highest 24 fim fluxes are safe disk detections.) For all objects with MIPS detections, 
the SEDs (given in flux units) are surprisingly flat: With three exceptions, which will be 
discussed below, the fluxes at 9-11 fim are consistent with the fluxes at 24 fim, and generally 
on t he level of a few m Jy. The MIR SED is sensitive to both disk structure and inclination 
(e.g. IWood et al.ll2002l ); since we can rule out high inclinations (i.e. close to edge-on disks) 



for all our targets based on the low visual extinction (see Sect. 13.11) . this points out that 
most brown dwarf disks in UpSco have similar structural characteristics, resulting in flat 
SEDs. 

The three exceptions are usco55, usdl61833, and usdl61939. These are the only objects 
with MIPS and IRS fluxes > 6mJy, and the only ones which show clear excess flux at 
24 fim, compared with 9 — 11 fim. Rising flux levels in the mid-infrared are a trade mark 
of st rongly flared disks, as show n in numerous comparisons o f observed and modeled SEDs 
(e.g. iKenyon fc Hartmannl 119871 ; iDullemond fc Dominikl 12004 ) . When the opening angle of 



the disk increases with its radius, the disk will intercept more light from the central object 
than a flat disk, leading to excessive mid-infrared emission. On the other hand, a flat disk 
where the dust has settled to the midplane causes a flat or gradually decreasing SED in the 
mid-infrared. Thus, the three objects with rising SEDs are probably the ones with the most 
flaring in our sample. Over all, the degree of flaring is low in UpSco brown dwarf disks, or 
conversely, most disks are likely affected by dust settling. 



at 



All three object s with rising SEDs are known to have companions ( see Table [J): u sco55 
; 18 AU (0.12") fcraus et al.ll2005h . usdl61833 at ~ 130 AU (0.9") jLuhmarJboOon . and 
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usdl61939 at ~ 25 AU (0.08" )Q faouv et al.ll2006bh . (All physical separations given in AU 



are projected values assuming a distance of 145 pc; because of the projection factor strictly 
speaken they should be considered as lower limits to the true separation.) The two closer 
companions are far from being resolved with MIPS and IRS (which have pixel sizes of 2.4 and 
1.8" respectively). The 240 AU companion of usdl618333 is not resolved by MIPS either, and 
it is likely to contaminate the IRS spectrum to some extent. Thus, our MIPS fluxes for these 
three objects include the flux from the companion. Photospheric fluxes from the companions, 
however, are not able to explain the 'flared' SED appearance, since the photospheres cannot 
contribute more than 0.5 mJy, less than 5% of our total 24 \xm fluxes for these three objects. 
Moreover, companions also contribute to near-infrared and IRS fluxes, thus cannot produce 
an increasing SED. Therefore, we attribute the excess emission at 24 fim to their disks. A 
discussion of possible influences of binarity on the disks follows in Sect. 15.41 

Comparable to the overall SED, the IRS spectra alone appear flat or decreasing in 
most cases. Only three objects - uscoll2, usdl60958, usdl61939 - show an indication of an 
emission feature (further investigated in Sect. 14.21 and Fig. [6]). Emission features around 
10 fim (often referred to as silicate features) are assumed to have its origin in an optically 
thin hot surface layer of small dust grains. The absence or weakness of this feature thus can 
be interpreted as lack of hot, small grains - possibly due to grain growth and subsequent 
dust settling to the (colder) disk midplane. Dust settling is confirmed for the majority of 
our objects by the general appearance of the SED (see above). Additional factors which can 
lead to dust settling and thus may weaken the silicate feature are a low gas-to-dust ratio or 
weak turbulence in the disk. We caution that disk inclination can also weaken the emission 
feature and thus might also play a role in some of our targets, although high inclinations are 
unlikely (see Sect. 13.11) . For a more detailed discussion of the minimal effects of inclination 
for our sample see also Fig. [5] and the assciated discussion in Sect. 14. 11 

Both findings, the uniformly flat SEDs and the paucity of strong dust emission features, 
demonstrate that UpSco brown dwarfs harbour disks that are clearly in an advanced stage 
of their evolution. At an age of 5Myr, flaring is not a major trait anymore. Instead, dust 
settling likely affects the majority of the disks, possible accompanied (and caused) by grain 
growth. We will come back to these findings in Sect. where we carry out comparisons 
with disk properties for younger and more massive objects. 



1 us dl61939 has additionally a wide companion which is, however, likely a background source ( Bouv et al 
2006bh . 
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3.3. Disk frequency 

In Fig. |5] we plot the K — 24 \xm colour vs. effective temperature for all targets. This plot 
will be used to derive an estimate for the disk frequency in our sample. Specifically, we want 
to find out for how many objects the measured MIPS flux (or upper limit) is consistent with 
photospheric emission. For this purpose, we calculate blackbody fluxes for the given effective 
temperature range at 2.2 and 24 /im, from which we obtain the photospheric K — 24fim 
colour. 

The solid line in Fig. [2] shows the result: Photospheric colours for our targets are 
between 0.9 and 1.4 mag, depending on effective temperature. Dashed lines show the colour 
assuming five times the photospheric emission at 24 fim. All objects with colours exceeding 
this level, have strong excess emission, and are thus likely to harbour a disk. Thirteen objects 
with MIPS detection are in this colour range, and are thus primary disk candidates. For one 
more object, the upper limit is too high to exclude excess due to a disk. For the remaining 
objects, the MIPS fluxes are fully consistent with pure photospheric emission. In total, the 
disk frequency for our sample is 37 ± 9% (13 out of 35, uncertainty corresponds to the la 
confidence interval based on binomial statistics), but might be as high as 40% (14 out of 35), 
if we count the upper limit above the threshold. 

This estimate is quite robust and does not depend on distance and radii of the ob- 
jects. Using a conservative limit (dashed line), we make sure that uncertainties in effective 
temperatures, which might be considerable (see above), and the inaccuracy of the black 
body approximation does not strongly affect our result. This conservative approach also 
prevents us from confusing disks with non-resolved, cool companions, which might enhance 
the mid-infrared flux. 

As can be seen in Fig. [2] there is a clear dichotomy in the mid-infrared colours of our 
targets: Down to T eff = 2600 K, objects with strong 24/im excess are separated by at least 
1.2 mag (i.e. a factor of three) from objects without excess. (This is a lower limit, because 
most objects without excess are in fact non-detections.) For the lowest mass objects, where 
the noise level increases, the difference between both populations is still at least a factor 
of two in colour. Thus, most of the objects with disks are clearly separated from objects 
without disks, leaving only a small number of faint objects with ambiguities. Thus, we 
can quite reliably conclude that, based on the 24 \im photometry, the disk frequency in our 
sample is 37%. This value, derived from an unbiased survey, is somewhat smaller, but still 
consistent with the disk frequencies obtained from smalle r samples of brown dwarfs in UpSco 



(50%, see IJayawardhana et al.ll2003at iBouy et al.ll2006al ). 



All 13 objects with MIPS detection and clear excess emission, i.e. objects with safe 



13 



disk detection, are marked in the last column in Table [TJ For clarity, they are also plotted 
with different symbols in Fig. [IJ Those are the primary candidates for the comparison with 
models, which will be described in Sect. HJ Their SEDs are shown in Fig. [3J 



4. Modeling the SEDs 

In this section, we will compare the observed SEDs with models to obtain more quan- 
titative constraints. Two different and independent approaches are used: a) a Monte Carlo 
radiative transfer model providing predictions for the MIR fluxes (Sect. 14.11 and b) linear 
combinations of mass absorption coefficients to fit the dust emission features around 10 /iin 
(Sect. 14. 2p . Both approaches give results consistent with the more qualitative discussion in 
Sect. 13.21 Brown dwarf disks in UpSco are ubiquitously affected by dust settling, possibly 
grain growth, and inner disk clearing. 



4.1. The overall SED 



We adopt the same modeling a pproach for the Up Sco Spitzer data as we did in our 
previous paper on brown dwarf disks (IScholz et al.ll2006l ). These SED models use our Monte 
Carlo radiatio n transfer code, ignore accretion, use N extGen model atmospheres for the 
input spectra (lAUard et al.ll200ll ; lHauschildt et al.lll999l ). and assume the disk surface density 
profile is given by a power law, S(-R) ~ R~ x . The disk density is parameterized by 



f =Po[ — I exp--[z/h(w)\' 



where w is the radial coordinate in the disk mid-plane and the scaleheight increases with 
radius, h = Hq (zu/R*) 13 . We vary the degree of flaring within the geometric disk models 
by adjusting the values of (3 and h . We assume two dust size distribut ions are present 



within the disk, a) small grains as with an interstellar-like size distribution (jKim et al.lll994l . 
typical grain size ~ 0.01 — 0.1 fim) and b) larger grains with sizes up to 1mm following 
a distribu tion we have used previously to model SEDs of disks around classical T Tauri 
stars (see IWood et al.l 120021 . for more details). The smaller, ISM-like grains have a larger 



scaleheight than a population of larger grains. This is a simple way to mimic dust coagulation 
and the settling of larger particles towards the disk midplane. In our parameterization we 
assign the mass M<j to the large grains and the mass /ismM* to the ISM-like grains. The dust 
mass is then (1 + /ism)^c2- F° r the conversion to the total disk mass, we assumed a gas to 
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dust ratio of 100. For all models w e assume that dust in regions close to the star is destroyed 
if temperatures rise above 1600 K (IDuschl et al.lll996l ). This condition provides a minimum 
inner dust radius of typically ~ 6i?*. Any remaining dust within this gap we assume to be 
optic ally thin and therefore we effectively have an opacity gap in the disk (ILada fc Adams 
19921). For further details of the model ingredients, assumptions, and dust properties see 
Scholz et all tood \. 



In our previous paper we discussed the many degeneracies in disk parameters that are 
inherent in SED modeling. In this paper we set the disk mass to be 4.5 x 1O _4 M . This 
disk mass is typical for the sources with 1.3 mm data that we modeled in our Taurus- Auriga 
sample. In the absence of mm data w e cannot constrain the disk mass, as the mid-IR SED 
is not very sensitive to disk mass (e.g. IWood et al.l 120021 ) . The lack of long wavelength data 



(sub-mm and mm) als o precludes us from d etermining the optical depths of the disks in our 
sample. As shown by I Walker et al.l (120041 ). low mass disks in hydrostatic equilibrium can 
provide acceptable fits for the mid-IR SEDs of brown dwarf disks. However, our models 
for Taurus brown dwarf disks showed that such low mass disks are not consistent with mm 
data. Therefore we believe that very low mass, optically t hin disks are unlike ly but mm 
data is required to test this further. As we demonstrated in lScholz et al.l (120061 ) our data is 
not sensitive to the outer disk radius, so we choose a disk radius of = 100 AU. For the 
object usdl60958, which has a close companion at ~ 12 AU and a disk (see Tabled]), we also 
tried a smaller disk radius of 5 AU. 

In Fig. H] we show the results of our SED fitting for the objects with disk; Table [2] 
contains the star and disk parameters used to model each source. Each SED shows the data, 
the input atmosphere model spectrum (dotted line), our model SED using the parameters 
fro m Table [2] (so l id line ) , and also the SED for a disk in hydrostatic equilibrium (dashed line, 



sec 



Scholz et al.l (120061 ) for a description of the hydrostatic disk models). For usdl60958, we 



plot additionally a model with disk radius of 5AU as dash-dotted line. All models are for 
a disk inclination angle i = 20°; none of our sources show the characteristic faint optical 
spectra indicative of highly inclined disks (see also Sect. 13. ip . 

As with our models for brown dwarf disks in Taurus, we immediately see that the highly 
flared hydrostatic disk models, which assumes that dust and gas are well mixed, produce in all 
cases larger IR excesses than observed. Thus, some degree of dust settling has to be assumed 
to match the observed IRS and MIPS fluxes. The models assuming two dust components 
without hydrostatic equilibrium agree much better with the observations, consistent with 
a picture of dust grain growth and the settling of larger grains towards the disk midplane, 
leaving a population of small grains suspended at large heights producing the silicate feature 
seen in a few cases (Sect. 13.21 and |4T2|) . The best fit is typically achieved with a small amount 
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of small grains or no small grains at all, i.e. when the dust mass is dominated by the larger 
grains. The SED of usdl 60958 can be modeled adequately assuming a small disk radius 
of 5AU, again demonstrating that the mid-infrared is not very sensitive to the global disk 
radius (see Sect. 15.41 for further discussion of the disk-binarity connection). 



As described in IScholz et al.l (120061 ). disk structure, inclination, and grain sizes are 
degenerate parameters in these models. To confirm that dust settling is required to match the 
observed SEDs, we explored these degeneracies in more detail in Fig. For one exemplarily 
case, the object uscoll2, we plot our best fit using the non-hydrostatic model (solid line 
for large and small grains, dashed line for small grains only) and hydrostatic, highly flared 
models (with small grains only) for different inclinations (dash-dotted lines). As can be seen 
in this figure, for low and moderate inclinations the hydrostatic models give clearly too much 
flux at 24/im, which turns out to be strongest constraint for the disk geometry. Conversely, 
we can exclude high disk inclinations (i.e. close to edge-on geometry) based on the low 
extinction seen in the near-infrared photometry. Thus, some extent of dust settling and/or 
grain growth is required to explain the observed SED, no matter what value is chosen for 
the inclination. Fig. [5] also demonstrates that non-hydrostatic models with predominantly 
large grains and with small grains only are both able to provide a reasonable match to the 
data, because the large-scale mid-infrared SED is not particularly sensitive to the grain size. 

The results from the model fits can be directly compared to our previous findings for 
younger brown dwarfs in Taurus, where we used the same models. In Taurus, however, we 
only analysed five out of 20 objects with detection at 1.3 mm; the sample is thus biased 
towards the most massive (and thus maybe less evolved) disks. For all five Taurus objects a 
small amount of ISM-like dust is required to fit the model SED. The same is found in UpSco: 
For five of our disks, the best fit is obtained when we add small grains. The values for the 
flaring parameters f3 are similar in both samples, but the scaling factors ho are much smaller 
(in absolute units) in all UpSco disks, leading to smaller scaleheights at a given radius. Thus, 
in our large sample in UpSco we do not see disks which are similarly strongly flared as the 
Taurus sources. This indicates a higher degree of dust settling in the UpSco objects when 
compared with younger brown dwarfs. 

As there is a gap in our data coverage in the 3 fim to 8 /im region, we cannot constrain 
the size of any inner disk holes in sources showing IRS and MIPS excesses. Sources without 
MIR excess, however, either have no disk or a disk with a inner hole, i.e. an optically 
thin region. To constrain the inner hole sizes, we applied the same models to five objects 
without excess emission (usco75, uscol04, uscol30, usdl555605, usdl62041). For four of 
these objects, the resulting model fits are shown in Fig. HI In all substantial inner 

hole is required to explain the lack of excess emission, the minimum hole sizes are 5AU for 
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uscol30, usdl55606, usdl62041, 7.5 AU for uscol04, and 20 AU for usco75. These estimates 
do not depend significantly on the choice of the grain size. 

Taken together, the results from the SED modeling provide strong evidence for dust 
settling in UpSco brown dwarf disks. Hydrostatic, highly flared models cannot reproduce 
the observed SEDs. In addition, we find that < 5 — 20 AU size (opacity) holes are required 
to fit the SEDs without disk excess. 



4.2. The silicate features 

The spectral region around 10 \xm is valuable for deriving the composition of the warm 
(T >200K) circumstellar dust. Indeed, it is here that the most common components of 
astronomical dust have spectral features. Using these features, d etailed analysis of the disk 



chemistry has bee n carried out previo usly for Herbig Ae/Be (e.g. iBouwman et al.ll200ll ) and 



4r 

T Tauri stars (e.g. iMeeus et all 120031 ). According to these studies, the main dust species in 
circumstellar disks are olivines, pyroxenes and silica, either in amorphous or crystalline state. 
A similar modeling approach, tailored to our faint targets, will be used in the following to 
analyse their IRS spectra. For a deeper discussion of the properties of astronomical dust 



and th e mass absorption coefficients k we used in this paper, we refer to Ivan Boekel et al. 
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Given that the flux levels around 10 /im are low for most of our targets, resulting in low 
signal-to-noise in the IRS spectra, we refrained from doing a detailed analysis of the dust 
composition, which requires high quality data. Instead, we tried to answer two questions 
using a more robust approach: a) Can we confirm the presence of silicate features in some 
of our targets (see Sect. 13.21) ? b) Is there evidence for grain growth from these features? 

Because changes in the dust composition only causes relatively subtle (for us unde- 
tectable) changes in the shape of the silicate feature, we simply adopted a dust composition, 
which contains only amorphous silicates, the most prominent dust species in many circum- 
stellar disks. To check for evidence of grain growth, we run the models with two different 
grain sizes, 0.1 and 1.5 /im. We represent the dust continuum by a blackbody with a temper- 
ature T cont , and for simplicity, assumed that all the grains emitting at 10 jum have the same 
temperature. Since we are modelling emission features, the region where the emission origi- 
nates must be optically thin. Therefore, the flux of the emission feature can be represented by 
a linear combination of the mass absorption coefficents: -F emisson (A) oc £>(T cont .) J2i a i K iWi 
where B{T cont ) is the Planck function with a temperature T cont ., and Kj(A) are the mass 
absorption coefficents and dj the multiplication factors for each individual species i (in our 
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case: only amorphous silicate). We used a reduced x 2 analysis to fit the spectra. For each 
object showing an indication of an emission feature, we obtain a best fit (lowest reduced x 2 ) 
for small and large grains, respectively. 

In the following, we summarise the results for our 13 sources with confirmed disk excess 
(see Sect. 13. 3p . For three objects, uscoll2, usdl60958, usdl61939, we detect a silicate 
feature, i.e. the spectrum can not be reproduced within the error bars only by a blackbody 
dust continuum. This is in line with the results from visual inspection of the spectra done 
in Sect. 13.21 For these three sources, we show the best fit for small and large grains in Fig. 
El As can be seen in the figure, both small and large grains give an adequate fit to the 
spectra. The x 2 values for the best fit do not differ significantly for small and large grains, 
respectively. Thus, even when using binned data, the flux levels and thus the signal-to-noise 
in our sample are too low to distinguish between different particle sizes from the IRS spectra. 

For the remaining 10 objects with disk the dust emission features around 10 [im are 
either absent or too weak to be confidently detected in our spectra. The MIR SEDs from 8 
to 13 fim are mostly flat and featureless. Since the dust emission has its origin in relatively 
small (< 10 fim), warm grains in the disk atmosphere, this implies that to some extent the 
dust has settled to the optically t hick midplane, a process wh ich is usually thought to be 



accompanied by grain growth (e.g. lDullemond &: Dominik!l2004l ). Disk inclination might also 



contribute in some cases to weaken the silicate feature, although this is not a major factor in 
our sample, because high inclinations (and thus close to edge-on disks) are excluded based 
on the low extinction (see Sect. I3.ip . 

In summary, we find silicate features for three out of 13 objects with disks. However, 
the flux levels and thus the signal-to-noise ratios are too low to give meaningful constraints 
on dust chemistry and particle sizes from the emission features. For the remaining 10 objects 
with disk we do not detect an emission feature. 



5. Discussion: Brown dwarf disks at 5 Myr 

Both the purely empirical discussion of the MIR SEDs in Sect. [3] and the results from 
the modeling (Sect. 14.11 and l4.2p provide a consistent picture of the disk properties in our 
sample. Given our large sample of coeval targets, they can be used to give constraints on 
evolutionary processes in brown dwarf disks. In the following, we will summarise the results 
from the previous sections and compare them with similar studies for younger objects and 
more massive stars (Sect. I5.ll and r5.2l) . Care has to be taken, though, when comparing results 
obtained at different wavelenghts, because longer wavelenghts probe the disk at larger radii. 
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This is further complicated by the fact that when changing central object masses the same 
wavelength probes different regions in the disks. Nevertheless, as it turns out, the main 
conclusions are basically unaffected by these restrictions. Furthermore, we will include a 
discussion of the accretion/disk connection (Sect. 15.31) and possible influences of binarity on 
disk evolution (Sect. 15.41) . 



5.1. Disk dissipation 



The main probe for disk dissipation timescales is the disk frequency. In Sect. 13.31 
we derive a disk frequency of 37 ± 9% for our sample of UpSco brown dwarfs, spectral 
types > M5, based on the 24 /im MIPS fluxes. This is consistent with values of 50% de- 
rived previously based on smaller samples in th e same region using shorter wavelength data 
( jJayawardhana et al.1 l2003al ; iBouy et al.l l2006al ). We can now compar e these results with 
similar studies for more massive stars. Recently, ICarpenter et al.1 (120061 ) published a Spitzer 
disk survey for a large sample of stars spanning a wide range in spectral types in UpSco, 
thus coeval to our targets. For stars with spectral types K0-M5 (corresponding to masses 
0.1-1 M , complementary to our targets), they report a disk frequency of 19 ± 4%§ derived 
from observations at 16 fim and a very similar value from observations at 8 fim. This is lower 
than our datapoint for brown dwarfs on a 1.8cr confidence level. 



Both samples, our brown dwarfs and the stars from lCarpenter et al.l (120061 ). are unbiased 
in the sense that the objects have been selected using criteria independent from their disk 
properties. The only difficulty when comparing the two studies is the different wavelength 
used to probe the presence of inner disks. However, all objects for which we see strong 
excess at 24 \im do also show significant excess already at IRS wavelengths, i.e. between 8 
and Yl\im. In addition, studies at wavelengths between 3 and 8 /im provide a simila r disk 
frequency (see above). Conversely, the disk fraction derived by lCarpenter et al.1 (120061 ) does 
not change significantly when going from 16 fim to 8 fim. Apparently, the disk fraction is not 
a strong function of wavelength between 8 and 24/im, regardless of object mass. Therefore, 
a comparison between the stellar disk frequency derived from 16 /im data and the substellar 
disk frequency obtained from 24 fim is valid. Moreover, both in the Carpenter et al. and 
in our sample the value for the disk frequency is fairly robust, because there is a clear gap 
in mid-infrared emission between objects with and without disks. Thus, this comparison 
suggests that the disk fraction for 5 Myr old brown dwarfs (spectral types >M5) is higher 



2 We re-calculated the uncertainty based on binomial statistics, to be consistent with the uncertainty 
estimate derived for our disk frequency. 
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(on a 1.8(T confidence level) than for coeval K0-M5 stars. 



The disk frequency (an d thus the average disk lifetime) in UpSco are therefore a function 
of stellar mass: As found by lCarpenter et all (120061 ). disks are less frequently observed around 
G/F stars compared with K0-M5 stars. Now we find that the disk frequency for objects with 
spectral types > M5 is even higher than for K0-M5. This indicates that the disk lifetime 
in UpSco, i.e. the age at which only very few objects (< 10%) retain an inner disk, is anti- 
correlated with the object mass: < 5Myr for G/F stars, but larger than 5Myr for brown 
dwarfs. 

To verify this result, it is useful to compare with IC348, a cluster with age of ~ 3 Myr, 



for wh ich disk frequencies have been published recently: Based on IRAC data, lLada et al. 
kood ) derive disk fractions of 11 ± 8% for stars earlier than K6, 47 ± 12% for K6 -M2 



stars, and 28 ± 5% for M2-M6 stars. Using a similar method, iLuhman et all (l2005bl ) has 
previously found a disk frequency of 42 ± 13% for the brown dwarfs (spectral types >M6) in 
this cluster. Taken together, these results give evidence for a mass- dependent disk lifetime, 
however, it may not simply be an increase towards lower masses as in UpSco. In pre-Spitzer 
studies, it was commonly f ound that brown dwarf disk lifetimes are not 'vastly different' 
(jJayawardhana et al.ll2003al ) from those of T Tau ri stars, i.e. for bot h object classes the disk 
frequency drops to values < 10% after 5-10 Myr (IHaisch et al.ll200ll ). While this conclusion 
still holds, Spitzer may now allow us t o determine disk li fetimes as a function o f mass more 
accurately. The results from IC3 48 (lLada et al.1 120061 ; ILuhman et all l2005bl ) and UpSco 
(this paper, ICarpenter et al.ll2006l ). although not entirely consistent, are first steps in this 
direction. Future studies in other star forming regions will hopefully help to disentangle the 
effects of mass, environment, and age on the disk dissipation timescale. 

Our data provide useful constraints for the duration of the inner disk clearing process. 
The brown dwarfs in UpSco clearly fall in two groups: About one third has significant MIR 
excess, whereas the MIR emission for the remaining objects is completely consistent with 
photospheric fluxes (see Sect. 13.31) . For the latter ones, we find that they either have no 
disks or disks with a large inner hole having a radius of at least 5-20 AU (see Sect. |4.ip . 
"Inner hole" in this context means an "opacity hole", i.e. the clearing from small (< 1mm) 
and thus observable dust grains; this does not necessarily imply that the region is completely 
devoid of material. 

This clear dichotomy and the lack of intermediate objects indicate that the dust clearing 
in the inner 5AU of the disks occurs rapidly. Assuming that our objects have an age spread 
of 2 Myr (see discussion in Sect. 13 . 1 j) and are evenly distributed in age, the time scale for 
this process is shorter than 2 Myr divided by the number of objects in our sample. Thus, a 
rough upper limit for the timescale of inner dust clearing within 5 AU is ~ 10 5 years. The 
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rapid clearing process is also supported by the fact that disk frequencies do not strongly 
vary when derived from observations in different wavelength regimes (i.e. from 3 to 24 /j,m, 
see above). 

The consequence is a pronounced scarcity of so-called transition disks, with lack of 
emission at A < 10 /im, but excess at 24 [im and thus a large hole in the inner disk. A 
few of these transition disks have bee n found recently and the y are believed to be 'caught 
in the act' of inner disk clearing (e.g. ICalvet et al.ll2002l . 120051 ). From our sample, we can 
constrain the frequency of transition disks among brown dwarfs in UpSco to be less than 8% 
(2a confidence limit). Since our wavelength coverage has a gap between 2.5 and 7.5 fim, this 
does not include objects with hole radii of < 1 AU. Therefore, our re sult is consistent with 
transition disk frequencies found for very low mass objects in IC34 8 (IMuzerolle et al.ll2006l ) 



as well as for stars in young clusters (e.g. ISicilia-Aguilar et al.l 120061 ). Based on the available 



data, we conclude that inner disk clearing timescales are not vastly different between stars 
and brown dwarfs. 

It should be emphasised that the inner disk clearing timescale is not comparable to the 
inner disk lifetime. In ner disk clearing appea rs to be a rapid process, which may be triggered 
by photoevaporation ([Alexander et al.ll2006l ) or formation of giant planets. It interrupts the 
relatively quiet viscous evolution of the disk, produces large inner opacity holes, and may 
initiate the transformation of T Tauri like accretion disks to debris disks. From all what we 
know today, this scenario does probably also apply in the substellar regime. 



5.2. Dust settling and grain growth 



Our results regarding dust evolution can be summarised as follows: Out of 13 definitive 
disk detections, 10 have a flat broad-band SED (measured in fluxes, within our uncertainties) 
from 8 to 24 fim, and only three show a weak flux increase at 24 \im. This is best explained 
by dust settling to the disk midplane, which is usually thought to be accompanied by grain 
growth. The radiative transfer modeling in Sect. 14. II confirms this result; for all objects with 
disk, some degree of dust settling is required to match the observed broad-band SED. 

For a model-independent evaluation of the geometric properties of the disks, we cal- 
culated the flux ratio F 2 4/s = F(24//m)/F(8//m) for the brown dwarf dis ks in UpSco 



(13 o bjects) and for five brown dwarfs in Taurus, recently investigated in IScholz et al. 
(120061 ). This flux ratio is a robust indicator for the flaring in the disk; it is expected to 



be around or below unity for flat disks and will exce ed unity significantly for strongly flared 



disks (e.g. IWood et al.l 120021 ; iMohanty et al.l l2Q04bl ) . From the observed SEDs, we obtain 
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-^24/8 = 0.5 ... 1.6 in UpSco and F 2 4/ 8 = 1.2 .. . 3.8 in Taurus. For three out of five objects 
in Taurus -F 2 4/ 8 exceeds the highest value in UpSco. Thus, our sample in UpSco lacks the 
highly flared disks typical for very young star forming regions, indicating a larger degree of 
dust settling. This confirms the results from comparing the outcomes of SED models for 
UpSco and Taurus discussed in Sect. 14.11 Again, it has to be emphasised that while our 
survey in UpSco is unbiased with respect to the disk properties, the small Taurus sample 
comprises only the objects with the most massive disks. 

While the dust settling can be probed reasonably well with our data, it is more difficult 
to evaluate the grain properties. For three objects in our sample, we find significant silicate 
emission features (see Sect. 14.2)) . indicating the presence of hot, small grains. Following 
Apai et al. J2005h (see their Fig. 2), we determined parameters for the shape (continuum- 
subtracted flux ratio F (11. 3 fim) /F (9. 8 fim)) and strength (peak flux above the continuum) 
for these three objects, using the continuum derived in Sect. 14.21 We obtain values of 0.5 to 
1.5 (±0.2) for the peak strength and 1.2 to 1.6 (±0.05) for the shape. Thes e numbers are no t 
strikingly different from w hat has bee n foun d for brown dwarfs in Chal ( lApai et al.l 120051 ). 
It has to be clarified that lApai et al.l (120051 ) have targeted objects known to be accreting 
and/or bright in the MIR. Thus, their sample is biased towards objects with strong MIR 
emission. Therefore, we basically compare the objects with the strongest silicate features in 
UpSco and Chal, and their dust properties (based on shape and strength of the emission 
feature) appear to be rather similar. This is based on small samples and should therefore be 
substantiated as soon as more data is available. To take into account the large number of 
objects without measurable silicate feature in UpSco, it would be particularly interesting to 
compare with a complete, unbiased brown dwarf IRS survey in a younger region like Chal 
or Taurus, which is not available yet. 

Neither the large-scale SED nor the silicate feature allow us to put much constraint on 
the coagulation of grains and their sizes. There are, however, two indirect arguments for the 
occurence of grain growth in UpSco disks: a) the evidence for dust settling, which is usually 
thought to happen as the grains coagulate and decouple from the gas, and b) the absence of 
strong silicate emission in the majority of the UpSco disks and thus the lack of hot, small 
dust grains. Based on these indirect arguments, the occurrence of grain growth is ubiquitous 
in our sample. 

The general picture of the dust evolution in brown dwarf disks by dust settling and 
grain growth is consistent with recent findings for stars. All these processes are now rou- 



tinely observed with Spitzer in T Tauri d isks (e.g. iMegeath et al.ll2005l ; iFurlan et al.l 12005 



Bouwman et al.ll2005l ; iForrest et al.l 120041 ) . It is not clear yet, however, to what extent the 
dust evolution is a function of object mass. This can now be probed in UpSco by comparing 
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the broad-band SEDs of stars and brown dwarfs. In Fig. [TJ we overplot the typical flux levels 
for K0-M5 stars with a nd without disks at 8 and 16 \xm as solid and dashed lines (taken from 



Carpenter et al.l (120061 )). Although the wavelength coverage for stars and brown dwarfs is 



not identical, hampering a reliable comparison, it is clear that the mid-infrared fluxes of the 
brown dwarfs with disks clearly fall in the range defined by the sample of stars. The spread 
in the mid-infrared SEDs appears to be smaller in the brown dwarf regime by 0.5-1 order 
of magnitude, maybe indicating a more uniform evolutionary timeline for very low mass ob- 
jects. But there is no clear evidence that the typical slope in the SED and thus the efficiency 
of the dust settli ng process i n the inner disk is a strong function of mass. This is contrary 
to the results of lLada et al.l (120061 ) in the somewhat younger cluster IC348, where a larger 
degree of dust settling is seen for late M-type objects. To what extent this discrepancy is a 
result of age effects or environmental differences has to be clarified in future work. 



5.3. Accretion signature from Ha emission 



We have compiled Ha 10% width measurements for our targets from our own spectra 
and the literature (see Table [p. This is the largest sample of brown dwarfs for which 
both MIR SEDs and Ha measurements are available. Ha emission is a useful indicator of 
disk accretion, and particulary th e 10% width has been found to be correlated with the 
accretion rate (INatta et al.l 120041 ) . Typically, a 10% width of ~ 300 kms -1 corresponds 
to accretion rates of ~ 10~ 10 M /year. A value of 200-270 kms -1 is usually adopted as 
a threshold betwe en accretors and non - accretors, corresponding to accretion rates below 
~ lO" n M /year Jwhite fe Basril 120031 : [javawardhana et al.l l2003bl ; iMohantv et~aD boosh . 
It is important to emphasise that this is mostly an empirical limit and in rare cases even a 
non-accretor might exceed the threshold because of fast rotation or strong magnetic activity. 

To probe the connection between disks and accretion, we plot the Ha 10% widths vs. 
the K — 24 \im colour already used in Sect. 13.31 (see Fig. [7]). The four quadrants in this 
figure correspond to the expected positions for objects without disks and without accretion 
(lower left), with disks but without accretion (lower right), with disks and accretion (upper 
right), and without disks but with accretion (upper left). The left half of this figure is easy 
to understand: There are no objects without mid-infrared excess and Ha 10% width clearly 
exceeding the 200 kms -1 limit. Thus, when an object has lost the disk or developed a inner 
opacity hole of at least 5 AU, the accretion also drops to unmeasurable rates. Thus, once the 
inner 5 AU of the disk is devoid of dust, the accretion rate drops below ~ 10~ n M /year. 

The right half of Fig. [7J is more puzzling. Surprisingly, only four out of 13 objects 
with clear evidence for the presence of an inner dust disk show a spectroscopic signature of 
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accretion. The remaining nine sources with inner disks are located well below the threshold 
between accretors and non-accretors, so their accretion rates are probably significantly lower 
than ~ 10~ n M /year. Thus, although a dusty disk is still present in the inner 5AU 
around these obje cts, accretion has ceased or at le ast reduced to a tr ickle. This confirms 
previous claims by I Jayawardhana et al.l (j2003al ) and lBouy et al.l (j2006al ) based on NIR/MIR 
observations of smaller samples. 

It is not clear to what extent the lack of accretion in these nine objects with disk is an 
indication for depletion of gas in the disk, particularly because two of these nine have the 
most flared disks in our sample. Flaring requires the presence of small grains (see Sect. 13.21 
and Sect. [HQ, and it is usually believed th at due to dust /gas coupling small dust grains are 



a good tracer for the presence of gas (e.g. lHaisch et al.l 1200 ll ). Thus, although accretion is 



not detectable in these seven objects, there is presumably still a (small?) amount of gas in 
their disks. 

A plausible explanation for the lacking disk/accretion connection might be variable ac- 
cretion. Spectroscopi c signatures of accretion a r e known to be varia ble in accreting objects 
(IJohns fc Basril Il995l ; lAlencar fc Batalhal 120021 ; iBouvier et al.l 12003k e.g.). Specifically for 
young brown dwarfs, we found recently that accretion rate changes by 0.5-1 order of _ mag- 



nitude and episodic phases of accretion and quiescence are not unusual (jScholz et al.l 12005 



Scholz &: Jayawardhana! 120061 ). Assuming that the phases of quiescence are dominant in the 
latest stages of accretion, maybe due to the onset of inner disk clearing, it might be that 
most objects with disks are observed in an epoch of low (and thus unmeasurable) accretion 
level, although they still have an inner disk. Multi-epoch spectra are required to constrain 
the degree of accretion variation in these objects. In summary, based on our large sample 
we find that the inner disks of brown dwarfs can persist longer than measurable rates of 
continuous accretion. 



5.4. Binarity and disk evolution 



Large parts of our sample have been observ ed with high spatial resolution, eithe r with 
HST or with AO imaging using NACO at VLT JBouv et al.ll2006bl ; IrCraus et alihoosh . The 
main goal of these surveys was to search for companions around young brown dwarfs; the 
observations are typically sensitive to substellar companions for separations between 5-10 AU 
and several hundred AU. The results are summarised in Tabled) Out of 29 objects that have 
been covered, 6 have been found to have a companion in the given separation range. These 
companions are unresolved in our Spitzer images and spectra. In the following, all separations 
are given as projected values, assuming a distance of 145 pc for UpSco. Strictly speaking, 
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these values should be considered to be lower limits to the physical separation. 



Two of the objects with companions in our sample, usco66 and uscol09, do not show 
any MIR excess, i.e. they have either lost their disk or the disk has a large inner hole. These 
are also the two objects with the smallest separations between the two binary components, 
10 and 5AU, respectively. The four remaining binary brown dwarfs (usco55, usdl60958, 
usdl61833, usdl61939) in the sample, with separations ranging from 12 AU to 130 AU, show 
clear signs for a dusty disk. 

The existence of companions inside 10 AU for objects without MIR excess is in line with 
the idea that these objects have either lost their disk or harbour a disk with a large inner 
hole, which would then be circumbinary. For the objects with inner disk and companion 
the situation is more complex. The presence of a companion is usually thought to truncate 
the individual disks of the components due to tidal fo rces, leaving two objects wit h small 
disks plus potentially an outer circumbinary disk (see lArtymowicz fc Lubowi 1 19941 ). which 
however would be undetectable with our MIR data. For solar-mass binaries with mass 
ratios close to one, the tidal trun cation radius is between 0.3 and 0.4 times the (physical) 
separation (lArmitage et al.lll999l ). For the closest binaries with disk - usdl60958, usco55, 
and usdl61939 - this might imply outer disk radii as low as 5-10 AU (assuming that projected 
separation is comparable to the physical separation). This scenario of disk truncation due 
to close companions is not inconsistent with our MIR SEDs: As we have shown in Sect. 14.11 
for the case of usdl60958 (which has a 12 AU companion), it is possible to have a small disk 
within the truncation radius that gives enough dust emission to explain the measured flux 
levels in the MIR (see Fig. |3j 7th panel). Thus, it appears plausible that these three objects 
feature small disks and companions. 

All three objects with increasing flux levels in the MIR (and thus flaring) have a com- 
panion, at projected separations of 18, 26, and 130 AU, as already pointed out in Sect. 13.21 
At face value, this might be a hint for enhanced flaring in the inner disk due to a compan- 
ion, e.g. by increased disk heating or turbulence. The literature, however, doe s not provide 



suppo rt for this scenario: Large-scale studies for T Tauri stars carried out by iFurlan et al. 



(120061 ) and iMcCabe et al.l (120061 ) find no correlation of MIR properties and the presence of 
companions. Thus, our finding of three objects with flared disks and companions might well 
be result of coincidence and should be treated with caution. 
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6. Summary 

We have carried out a Spitzer survey of 35 brown dwarfs in the Upper Scorpius star 
forming region with an age of ~ 5 Myr. Combined with literature data from spectroscopy 
and high-resolution imaging, the brown dwarfs in UpSco now comprise a large sample of 
substellar objects, for which information about disks, accretion, and multiplicity is available. 
In the following, we list our most important results: 



Out of 35 objects, 13 show clear excess emission at 24 fim, indicating emission from a 
dusty circum-sub-stellar disk. This corresponds to a disk frequency of 37 ± 9%. The 
disk frequency is on a 1.8a c onfidence level higher i n our sample than for K0-M5 stars 



in the same region (19 ± 4% ICarpenter et al.ll2006l ). suggesting a mass-dependency in 
the disk lifetime. 

2. The objects without emission either do not have a disk or they harbour a disk with 
large inner holes (radii > 5AU). The clear distinction between objects with/without 
disks and the lack of objects in transition between the two subsamples is a clear sign 
of rapid inner disk clearing, probably on timescales of < 10 5 years. 

3. The midinfrared SEDs of the disks appear uniform: Most of them are flat between 8 
and 24 fim with flux levels of a few mJy Only three objects show slighly enhanced 
flux levels at 24 /im, indicative of flaring. Hydrostatic disk models where dust and gas 
are well-mixed are not able to fit the data. However, assuming a two component dust 
composition where large grains have smaller scaleheights than ISM-like small grains 
provides good agreement between model and observations. This is well-explained with 
a scenario in which grain growth leads to larger grains that settle to the disk midplane. 

4. In most disks, there is no clear indication for dust emission features around 10 /im, 
arguing for a lack of hot, small grains. This confirms that the disks are affected by 
dust settling and grain growth. 

5. Brown dwarf disks in UpSco (5 Myr) show a larger degree of dust settling than in 
the younger Taurus star forming region. In UpSco itself, we do not see a significant 
difference between stars and brown dwarfs in the broad-band SEDs, indicating that 
dust settling is not a strong function of object mass. 

6. Taken together, our survey finds disks at an advanced stage: A large fraction of the 
objects have already cleared out inner parts of the disk. Evolutionary processes like 
dust settling and grain growth are ubiquitous in brown dwarf disks at 5 Myr. Since 
these processes are prerequisites for the formation of larger bodies in the disk, our 
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results suggest possible planetesimal growth and planet formation in brown dwarf 
disks. 

7. Out of 13 objects with disks, only four appear to accrete at more than ~ 1CT 11 M /year. 
Thus, dusty disks can persist although continuous accretion has stopped or at least 
reduced to unmeasurable levels. 

8. Four objects with disks also have a companion; in three cases the (projected) separa- 
tions are smaller than 30 AU. This can be explained with the existence of a small disk 
(radius 5-10 AU) around the primary, truncated by the companion. 

This work is based on observations made with the Spitzer Space Telescope in the frame- 
work of the GO program #20435 (PI: Ray Jayawardhana) and the GTO program #248 (PI: 
Subhanjoy Mohanty). We thank Subhanjoy Mohanty for unbureaucratic data exchance and 
stimulating discussions. The help of Duy Cuong Nguyen, who was the observer at Magellan, 
is gratefully acknowledged. Fruitful discussions with the star formation group in Toronto, 
in particular with Yanqin Wu, have helped to improve the paper. We thank the referee for 
a thorough and constructive report. This project makes use of data products from the Two 
Micron All Sky Survey, and was supported by an NSERC grant and University of Toronto 
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Table 1. Object parameters and mid- infrared fluxes. 
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0.68 


0.29 


< 0.4 


usdl60958 


16 


09 


58.5 


-23 45 18.6 


M6.5 


2900 


363 


12 


4.32 


4.23 


4.03 


0.24 


3.73±0.37 b 


usdl61005 


16 


10 


05.4 


-19 19 36.0 


M7 


2850 


228 


s 


2.81 


2.23 


1.99 


0.51 


3.04±0.30 b 


usdl61006 


16 


10 


06.0 


-21 27 44.6 


M8.5 


2500 


110 c 


NaN 


0.91 


1.52 


1.62 


0.54 


0.78±0.20 b 


usdl61007 


16 


10 


07.5 


-18 10 56.4 


M6 


3000 


89 


S 


2.27 


1.79 


1.01 


0.60 


< 0.4 


usdl61103 


16 


11 


03.6 


-24 26 42.9 


M9 


2400 


90 c 


NaN 


1.61 


1.61 


1.31 


0.34 


0.99±0.20 b 


usdl61452 


16 


11 


52.6 


-20 17 13.2 


M9 


2400 


79 c 


NaN 


0.15 


0.15 


-0.02 


0.11 


< 0.4 


usdl61632 


16 


16 


32.2 


-22 05 20.2 


M6 


3000 


93 


S 


1.83 


1.68 


0.95 


0.52 


0.35 ±0.20 


usdl61833 


16 


18 


33.2 


-25 17 50.4 


M6 


3000 


72 


130 


5.35 


4.37 


4.16 


0.33 


7.47±0.75 b 


usdl61840 


16 


18 


40.8 


-22 09 48.1 


M7 


2850 


90 


S 


0.70 


0.29 


0.51 


0.51 


< 0.4 


usdl61903 


16 


19 


03.4 


-23 44 08.8 


M6.5 


2900 


94 


s 


0.76 


0.86 


0.78 


0.28 


< 0.4 


usdl61916 


16 


19 


16.5 


-23 47 22.9 


M8 


2600 


93 c 


NaN 


1.06 


0.88 


0.65 


0.45 


0.5-l b 


usdl61926 


16 


19 


26.4 


-24 12 44.5 


M6 


3000 


161 


S 


0.80 


0.84 


0.16 


0.27 


< 0.4 


usdl61929 


16 


19 


29.9 


-24 40 47.1 


M8 


2600 


135 c 


S 


0.49 


0.92 


0.41 


0.29 


< 0.4 


usdl61939 


16 


19 


39.8 


-21 45 35.1 


M7 


2750 


260 


26 


7.60 


8.53 


8.78 


0.52 


10.7 ± l.l b 


usdl62041 


16 


20 


41.5 


-24 25 49.0 


M7.5 


2650 


81 c 


S 


0.55 


0.49 


0.44 


0.41 


< 0.4 



a Upper limits are based on 2cr flux upper limits. Errors correspond to la flux uncertainties. 
b Safe disk detections according to the analysis in Sect. 13.31 
c From Ijay awardhan a et al. | (|2002ll and lMohanty et al. | ll2005h . 

d Binarity: If a companion has been found, the projected separation in AU is given, assuming a distance of 145 pc. A 'S' ('single') 
indicates that the obj ect does not have a companion at se parations > 5AU. Objects with 'NaN' have not been observed yet with high 
spatial resolution fsee lBouv et al. 2006b; Krau set al.ll20o'il) . 

°Also called usdl60019, see lMartm et al.1 l l2004h . 
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Table 2. Parameters from the radiative transfer modeling of the overall SED (see Sect. 

ED 



Full name 




R* 




flout 


c 


hf m 


/ism 


P 




(K) 


(Re) 


(Mq) 


(AU) 


(R*) 


(R*) 






usco55 


2800 


0.65 


0.26 


100 


0.015 




0.0 


1.15 


uscoll2 


2850 


0.40 


0.10 


100 


0.01 


0.02 


0.03 


1.1 


uscol28 


2600 


0.30 


0.01 


100 


0.0375 




0.0 


1.1 


usdl55556 


2900 


0.40 


0.075 


100 


0.06 




0.0 


1.0 


usdl55601 


2900 


0.35 


0.075 


100 


0.03 




0.0 


1.1 


usdl60603 


2900 


0.35 


0.03 


100 


0.05 




0.0 


1.0 


usdl60958 


2900 


0.55 


0.075 


100 


0.01 




0.0 


1.1 




2900 


0.55 


0.075 


5 


0.007 




0.0 


1.1 


usdl61005 


2850 


0.33 


0.04 


100 


0.025 


0.03 


0.03 


1.1 


usdl61006 


2500 


0.26 


0.02 


100 


0.05 


0.06 


0.03 


1.0 


usdl61103 


2500 


0.26 


0.02 


100 


0.05 


0.06 


0.03 


1.0 


usdl61833 


3000 


0.65 


0.10 


100 


0.0055 




0.0 


1.2 


usdl61916 


2600 


0.25 


0.03 


100 


0.015 


0.025 


0.03 


1.1 


usdl61939 


2750 


0.50 


0.04 


100 


0.03 


0.04 


0.03 


1.15 
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Fig. 1. — Broadband spectral energy distributions for all our sources, normalized to the J 
band flux. Objects with excess at 24 \xm (see Fig. [2]) are shown as filled squares, objects 
without excess as plus symbols. The arrow at 24 fim gives a typical upper limit for objects 
without MIPS detection. Dashed lines show average flux levels for K0- M5 stars without 
disks in UpSco, solid lines the typical range for stars with disks (taken from lCarpenter et al. 

hood }). 
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Fig. 2. — K-band minus 24/zm colours for UpSco brown dwarfs calculated from MIPS fluxes 
vs. effective temperature. Upper limits correspond to 2a upper limits at 24/xm. The solid 
line show the photospheric colours estimated by assuming blackbody radiation; the dashed 
lines show colours assuming five times the photospheric emission at 24/zm and thus significant 
excess (see Sect. 13 .3 j) . 
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Fig. 3. — Spectral energy distribution for objects with MIR excess emission in comparison 
with models (dotted lines: photospheric flux, dashed lines: photosphere + hydrostatic disk, 
solid lines: photosphere + two-component disk, see Sect. 14.11) . For usdl60958, a model with 
disk radius of 5 AU is shown as dash-dotted line (7th panel). The parameters for the models 
shown in solid lines are listed in Table EJ 
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Fig. 4. — Spectral energy distribution for objects without 24 fim MIR excess emission in com- 
parison with models (dotted lines: photospheric flux, solid lines: photosphere + disk). The 
SEDs can be explained by photospheric emission plus disks with large inner holes (holesizes 
> 5AU; see Sect. SJJ. 
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Fig. 5. — Detailed SED modeling for the object uscoll2: The solid, thick line shows our 
best (non-hydrostatic) model already plotted in Fig. [3] (2nd panel), using small and large 
dust grains. The dashed line is a best fit using the same model approach, but with small 
grains only. Dash-dotted lines show hydrostatic models (with small grains only) viewed at 
ten different inclinations evenly spaced in cos(i) from i=0 to i=90. The photosphere is again 
shown as dotted line. 
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Fig. 6. — IRS spectra for three objects with detection of dust emission features, together 
with uncertainties for the individual datapoints. Additional calibration uncertainties are 
~ 10%. The grey lines show a model fit to the feature using either small (left panel) or large 
grains (right panel), see Sect. 14.21 for details of the modeling). As can be seen, both grain 
species provide a good fit to the data. The black solid lines show the continuum from the 
fit. 
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Fig. 7. — Ha 10% width vs. K-band minus 24/im colours for UpSco brown dwarfs. Upper 
limits correspond to la upper limits at 24/im. The dashed horizontal line separates accretors 
from non-accretors, while the dashed vertical line separates objects with infrared excess from 
those without (see Sect. 15.31) . 



